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ABSTRACT 

We combine two Gaussianization techniques - Wavelet Non-Linear Wiener Filter 
(WNLWF) and density reconstruction - to quantify the recovery of Fisher information 
that is lost in the gravitational collapse. We compute a displacement fields, in analogy 
with the Zel'dovich approximation, and apply a Wavelet Non-Linear Wiener Filter 
that decomposes the reconstructed density fields into a Gaussian and a non-Gaussian 
component. From a series of 200 realizations of N-hody simulations, we compute the 
recovery performance for density fields obtained with both dark matter particles and 
haloes. We find that the height of the Fisher information trans-linear plateau is in- 
creased by more than an order of magnitude at k > l.Oft-Mpc -1 for particles, whereas 
either technique alone offers an individual recovery boost of only a factor of three to 
five. We conclude that these two techniques work in a symbiosis, as their combined 
performance is stronger than the sum of their individual contribution. When applied 
to the halo catalogues, we find that the reconstruction has only a weak effect on the re- 
covery of Fisher Information, while the non-linear wavelet filter boosts the information 
by about a factor of five. We also observe that non-Gaussian Poisson noise saturates 
the Fisher information, and that shot noise subtracted measurements exhibit a milder 
information recovery. 

Key words: cosmology: theory — dark matter — large scale structure of universe — 
methods: statistical 



1 INTRODUCTION 

Understanding the nature of dark energy has been iden- 
tified internatio nally as one of the main goal of mod- 
ern cosmology (|Albrecht et alj [2006), and many dedi- 
cated experime nts attempt to constrain its equation o f 
state: LSS tFI (ILSST Science Collaborations et all [20091) . 
EUCLIE0 dBeaulieu et al l l2010h . JDElVl 3 ! jGehrehl hoioh . 
CHIMrjT [Peterson et aTT 120061 ) , SKAQ (ISchilizzil 120071 ; 
iDewdnev et all l2009l h BOSSFI (|Schlegel et all |2009| 1 and 
Pan-STARR^j. One of the favored technique involves a de- 
tection of the Baryonic Acoustic Oscillations (BAO) sig- 



nal dSeo fc Eisensteinl 12001 120051 : lEisenstein et all 120061 : 
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tection of the BAO wiggles in the matter power spectrum, 
which act as a standard ruler and allow one to map the 
cosmic expansion. 

With the new and upcoming generation of dark en- 
ergy experiments, the precision at which we will be able 
to measure the cosmological parameters is expected to drop 
at the sub-percent level, therefore it is essential to under- 
stand and suppress every sources of systematic uncertainty. 
In a BAO analysis, one of the main challenge is to extract 
an optimal and unbiased observed power spectrum, along 
with its uncertainty; the latter propagates di rectly on the 
dark energy parame ters with Fisher matrices (|Fisher1 [19351 ; 
iTegmark et al.l ll997). This task is difficult for a number of 
reasons. 

For instance, the scales that are relevant for the anal- 
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yses sit at the transition between the linear and the non- 
linear regime, at least for the redshift at which current 
galaxy surveys are sensitive, hence the underlying uncer- 
tainty on the matter power spectrum is affected by the 
non-linear dynamics. Th ese effectively coup les the phases of 
different Fourier modes l) Zhang et al ] l2003h and the Gaus- 
sian description of the den s ity fields has been observed to 
fail dMeiksin fc Whitdl 19991; iRimes fc Hamiltor]|2005l . 120061 ; 
iNevrinck et al.ll2006l : lNevrinck fc Szapudill2007T ). For an esti- 
mate of the BAO dilation scale to be robust, one must there- 
fore include in the analysis the full non-lin ear covariance of 
the po wer spectrum. Although results fromlTakahashi ~et al.l 
(2011) seem to suggest that non-Gaussianities had no real 
effect on the final results, it was recently shown that this 
was only true if the original power spectrum was measured 
in an unbiased an d optimal way, which is rarely the case 
l|Ngan et al.ll201lh . Otherwise, the discrepancy on the con- 
straining power is at the percent level. One of the way to 
reduce the impact of the non-linear dynamics is to trans- 
form the observed field into something that is more linear. 
Over the last few years, many "Gaussianization" techniques 
have been developed, which all attempt to undo the phase 
coupling between Fourier modes. The number of degrees of 
freedom - i.e. uncoupled phases - can be simply quanti- 
fied by the Fisher information, and recovering parts of this 
erased information can lead to improvements by factors of 
a few on cosmological parameters. 

For example, a de nsity r econstructi on al- 
gorithm ilEisenstein et al.1 120071 ; iNoh et all 120091 ; 



iPadmanabhan et al.l 2009T ). based on the Zel'dovich 



approximation, has been shown to reduce b y a factor of two 
the constraints on the BAO dilation scale ( Eisenstein et all 



120071 ; iNgan et"al] l201lh . This technique was applied on 
the SDSS data recently (|Padmanabhan et alj I2012T ) to 
improve the BAO detection, with small modifications to 
the algorithm such as to correct for the survey selection 
function and redshift space distortions. As discussed 
therein, an important issue is that two main mechanism are 
reducing our ability to measure the BAO ring accurately: 
1) a large coherent ~ 50 Mpc infall of the galaxies on 
to overdensities, which tends to widen the BAO peak, 
and 2) local non-linear effects, including non-linear noise, 
which also erase the smallest BAO wiggles. Reconstruction 
addresses the first of these mechanisms, and its is important 
to know whether something can be done about the second, 
after reconstruction has been applied. 

Wavelet Non-linear Wiener Filters (hereafter WNLWF, 
or just wavelet filter) were used to decompose dark mat- 
ter density fields dZhang et al.ll201ll ) and weak gravitational 
lensing K-fields |Yu et al. 2012T ) into Gaussian and non- 
Gaussian parts, such as to condense in the latter most of 
the collapsed structure; the Gaussian part was then shown 
to contain several times more Fisher information than the 
original field. Ot her method s include log-normal t ransforms 
<|Seo et alJkoill ). Cox-Box (Ijoachimi et al.ll201ll). running 
iV-body simulation backwards ( Goldberg fc Spergelll200oh . 
or dir ect Gaussianiza tion of the one-point probability func- 
tion |Yu et al.l l201ll ). just to name a few. This technique 
seems perfectly suited to address the issue of non-Gaussian 
noise described above. 

Our focus, in this paper, is to discuss how two of these 
techniques can be used in conjunction to maximize the re- 



covery of Fisher information. Not all combinations of Gaus- 
sianizat ion techniques a re winning, however, ft was recently 
shown i|Yu et all [2012T ) that WNLWF and log-transforms 
are not combining in an advantageous way. On one hand, 
if the log-transform is applied onto a Gaussianized field, the 
prior on the density field is no longer valid, and the log- 
transform maps the density into something even less Gaus- 
sian. On the other hand, it was shown that the log-transform 
is less effective than WNLWF alone at recovering Fisher in- 
formation, at least on small scales. Applying the filter after 
the log-transform does not improve the situation, since the 
Gaussian/non-Gaussian decomposition is less effective. In 
other words, the Fisher information that the log-transform 
could not extract is not recovered by WNLWF, and we are 
better off with the WNLWF alone. 

It seems, however, that this unfortunate interaction is 
not a constant across all combinations. In this paper, we dis- 
cuss how non-linear Wiener filters, constructed in Wavelet 
space, can improve the results of a density reconstruction 
algorithm, which takes the density back in time using linear 
perturbation theory. Our first result is that these two tech- 
niques work well together, in the sense that the final Fisher 
information recovery is larger than the two techniques stand 
alone. We first obtain these results with particle catalogues 
extracted from A r -body simulations, and extend our tech- 
niques to halo catalogues, which provide a sampling of the 
underlying matter field that is much closer to actual obser- 
vations. 

The structure of the paper is as follows: in Section [3] 
we briefly review the theoretical background of the density 
reconstruction and WNLWF, and review how we extract 
the density power spectra, their covariance matrices, and 
the Fisher information. We discuss our results in Section [3] 
and conclude in Section fj] 



2 THEORETICAL BACKGROUND 
2.1 Numerical Simulations 

Our sample of 200 A r -body simulations are gen erated with 
CUB EP3M, an enhanced version of PMFAST (|Merz et all 
2005) that solves Poisson equation with sub-grid resolution, 
thanks to the p 3 m calculation. Each run evolves 512 3 par- 
ticles on a 1024 3 grid, and is computed on a single IBM 
node of the Tightly Coupled System on SciNet (|Loken et alj 
120101 ) with n M = 0.279, Q A = 0.721, a 8 = 0.815, n s = 0.96 
and h — 0.701. We assumed a flat universe, and started 
the simulations at Zi — 50. Each simulation has a side 
of 322.36ft _1 Mpc, and we output the particle catalogue at 
z = 0.054. We search for haloes with a spherical over-density 
algorithm (|Cole fc Lacevlll996l ) executed at run time, which 
sorts the local grid density maxima in descending order of 
peak height, then loops over the cells surrounding the peak 
center and accumulates the mass until the integrated den- 
sity drops under the collapse threshold of 178, and finally 
empties the contributing grid cells before continuing with 
the next candidate, ensuring that each particle contributes 
to a single halo. Halo candidates must consist of at least 
one hundred particles, ensuring the haloes are large and col- 
lapsed objects. The center-of-mass of each halo is calculated 
and used as position, as opposed to its peak location, even 
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though both quantities differ by a small amount. We men- 
tion here that algorithms of this kind have the unfortunate 
consequence to create an exclusion region around each halo 
candidate, thus effectively reducing the resolution at which 
the halo distributions are reliable. Each field contains about 
88, 000 haloes, for a density of 2.6 x 10~ 3 /i 3 Mpc~ 3 . For com- 
parison, this is abo ut eight times large r than the density of 
the BOSS density l|Schlegel et alj|2009h . 

2.2 Density reconstruction algorithm 

We use a density reconstruction algorithm that is based 
on the linear the ory prediction, first found by Zel'dovich 
(|ZerDovichlll97oT ). that couples the density field <5(q, to) to 
the displacement field s(q) via 



f(q,to) = - V-s(q) 



(1) 



In the above expression, q is the grid, or Lagrangian, co- 
ordinate, and the displacement field is obtained in Fourier 
space as 



bQl) = ~8(k,to)F{k) 



(2) 



where F(k) = exp[— (kR) 2 /2] is a smoothing function sup- 
pressing features smaller than R = 10h~ Mpc. Particles at 
Eulerian coordinate x are displaced from their grid positions 
following 



x(t Q ) = o(t)[q + D(to)s(q)] 



(3) 



with D(t) the linear growth factor. These calculations are 
commonly used for the generation of initial conditions in 
iV-body simulations, and are accurate as long as the small- 
est scales probed are still on linear regime at the starting 
redshift. In the case where the particles - or haloes - to be 
displaced are not at to, one must subtract from the result 
the di splacement field from the grid location fsee lNoh et al.l 
(2009) for a detailed explanation of this technique). 

2.3 Wavelet non-linear Wiener filter (WNLWF) 

In this subsection we b riefly review the WNLWF algorithm , 
and direct the reader to lZhang et~ai1 (|201lMYu et all (|2012l ) 
for more details. 

W e consider in this paper the Daubechies-4 (|Daubechiesl 
Il992h discrete wavelet transform (DWT), which contains 
certain families of scaling functions (f> and difference func- 
tions (or wavelet functions) ip. The density fields are ex- 
panded into combinations of these orthogonal bases, and 
weighted by scaling function coefficients (SFCs) and wavelet 
function coefficients (WFCs). In our WNLWF algorithm, we 
deal with only WFCs, each of which characterizes the am- 
plitude of the perturbation on a certain wavelength and at 
a certain locations. 

In the three dimensional case, the properties of each 
perturbation depend on three scale indices (31,32,33) - con- 
trolling the scales of the wavelet Daubechies-4 functions 
- and three location indexes (h,fo,h) - controlling their 
translations. Specifically, on a given dimension, the grid scale 
corresponding to a specified dilation is L/2 3 (L — 1024 
in our case), and the spatial location is determined by 
IL/2 1 ' < x < (l + l)/2 3 . After the wavelet transform, all SFCs 
and WFCs are stored in a 3-dimensional field, preserving 



the gr id resolution (see iFang fc Thews] l| 19981 '); Ipress et all 
(1992) for more details). 

Our non-linear Wiener filter (NLWF) construction 
strategy relies on the fact that in wavelet space, the non- 
Gaussianities are clearly characterized in the probability dis- 
tribution function (PDF) of the WFCs e 31 J2i j 3 -i y ,i 2 ,i 3 . We 
thus construct our filter by splitting the wavelet transform 
of the original density, which we label D, into a Gaussian (G) 
and a non-Gaussian (JV) contribution. Since wavelet trans- 
forms are linear operations, this Gaussian/non-Gaussian de- 
composition happens also in real space when we wavelet 
transform back the contributions. Namely, we can write in 
wavelet space and real space respectively, 



D — G + N and d = d G + d N G 



(4) 



where the original density (d) is expressed as the sum over 
a Gaussian contribution (do) and a non-Gaussianized con- 
tribution (dNG). Our goal is thus to design a filter that con- 
centrates most of the collapsed structure in dNG, and thus 
produces dG that are closer to linear theory. 

The NLWF acts on individual wavelet modes, which are 
defined as combinations (not permutations) of all WFCs 
having the same three scale indices (ji, 32, 33)- For each 
wavelet mode, the NLWF is determined completely by the 
PDF f(x) of the corresponding WFCs, which is constructed 
by looping over the other three indices (h, l2,h)- We then fit 
this P DF with the analytical function presented in lYu et al.l 
(2012) (equation (15)): 



/pdf(i) 
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r(i as 2 



-0 



(5) 



and extract the two parameters a and s. These are actually 
dependent on the second and the fourth central moment of 
the PDF f(x) TO2 and m.4, hence we measure the moments 
first, then extract a and s via: 
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We then loop back over all spatial indices {l\,l2,lz) of this 
wavelet mode and decompose each WFC into two compo- 
nents: 



(ln/)» / x 



wng(i) 



(m /)'(*) 



1 + 



(7) 



(8) 



which are functions of x only. Note that the final filter func- 
tion depends only on s, which corresponds to the full width 
at half maximum of the Gaussian NLWF function iug. It 
characterizes the extent of the departure from a Gaussian 
PDF: the greater the s, the smaller departure is from Gaus- 
sian statistics. In comparison, a parameterizes the central 
deviation of the PDF: cr ce ntrai = a~ 2 . The same decompo- 
sition is performed on the reconstructed density fields and 
on those obtained from the halo catalogues. In this paper, 
we do not make use of the information contained in the 
non-Gaussian component and simply discard it, although it 
serves as a powerful probe of small scale structures and could 
help identifying haloes in a (Gaussian) noisy environment. 
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2.4 Information recovery 

The calculation of uncertainty about dark energy cosmolog- 
ical parameters is based on a propagation of the uncertainty 
about the matter power spectrum. In this process, the num- 
ber of degrees of freedom - i.e. the Fisher information - 
contained in the field is directly related to the constraining 
power. In this section, we review how the Fisher informa- 
tion about the amplitude of the matter power spectrum is 
calculated from simulated dark matter particles and haloes. 

The power spectrum P(k) of a density contrast <5(x) is 
calculated in a standard way: 



P(k) = <P(k)) = (\5(k)\ 2 ) 



(9) 



where the angle brackets refer to an average over our 200 
different realizations and over the solid angle. 

The uncertainty about the power spectrum is estimated 
from a covariance matrix C is defined as 



C(k,k') 



N 



1 N 

~ (P(k))][Pi(k') ~ (Pffl, (10) 



where N is the number of realizations and (P(k)) is the mean 
angular power spectrum over all realizations. The cross- 
correlation coefficient matrix is somehow more convenient 
to plot since it has higher contrasts, and is defined as: 



p{k,k') 



C(k,k') 



y/C(k,k)C(k',k') 



(11) 



The diagonal is normalized to one, and each element rep- 
resents the degree of correlation between the scales (k,k'). 
In the Gaussian approximation, the density is completely 
characterized by the power spectrum. Namely, 



C G {k,k') 



2P 2 (fc) 
N(k) ' 



(12) 



Consequently, pa is identical to the identity matrix. As ex- 
pected from the theory of structure formation, the non- 
linear collapse of the matter density field tends to couple 
the Fourier modes, which are otherwise independent, start- 
ing from the smallest scales and progressing towards larger 
scales with time. This coupling is responsible for highly cor- 
related region in p, and can be understood in terms of higher- 
order corrections, including the bispectrum, the trispectrum, 
etc. 

The Fisher information measures the number of inde- 
pendent Fourier modes in a density field up to a resolution 
scale fc max . We see from equation (|12p that dividing the co- 
variance by the (square of the) power spectrum is propor- 
tional to the number of independent measurements N(k). 
Therefore, the normalized covariance is defined as: 



C n < u rn (fc,fc') 



C(k,k') 
P(k)P(k') 



(13) 



Then, the number of degrees of freedom up to a scale fc max 
is otained by inverting the corresponding sub-sample of the 
normalized matrix, then summing over all the elements: 



/(fcn 



(14) 



The inversion of the covariance matrix involved in the cal- 
culation of the Fisher information amplifies the noise, hence 



such measurements typically requires a very strong conver- 
gence on the forward matrix. This can otherwise lead to 
biases of a fe w percent on deriv ed quantities like the BAO 
dilation scale (|Ngan et al.ll201lh . Generally, a forward ma- 
trix that is closer to diagonal contains more Fisher informa- 
tion; the theoretical maximum corresponds to the Gaussian 
case, where all modes are independent. Any non-vanishing 
off-diagonal element reduces the information. 



3 RESULTS 

In this section, we describe and quantify our ability at re- 
covering Fisher information with our two Gaussianization 
techniques, for density fields measured with simulated par- 
ticles and haloes. We recall that the later is much closer 
to actual observations, since galaxies trace highly collapsed 
structures. 



3.1 Density fields 

To illustrate the effect of different Gaussianization tech- 
niques, we present in Fig.[T]the projections through a thick- 
ness of 50 cells of a given realization after density recon- 
struction alone, after WNLWF alone, and with both tech- 
niques applied. We observe that the reconstruction reduces 
the size of each halo, as expected from this algorithm: par- 
ticles attempt to travel out of the gravitational potential. 
WNLWF has a slightly different visual effect on the den- 
sity: it removes most of the smallest structure p erturbations , 
leaving behind the larger ones. As discussed in [Penl i|l999l ): 
IZhang et all (|201ll ). the geometry of the Cartesian wavelet 
leaves behind a grid patterns, which only affects the small- 
est scales of the power spectrum and has no impact on the 
scales we are interested in. The combination of both tech- 
niques is presented in the middle right panel, and visually 
presents the least amount of collapsed structures. We also 
show the non- Gaussian part of the wavelet filter with and 
without reconstruction in the bottom panels. The largest 
peaks and sharpest structures are indeed filtered out. 

Fig. [2] shows the power spectrum of the dark matter 
particles and haloes before and after the Gaussianization 
techniques. We first observe that the measurement form the 
original particle field agrees at the few percent l evel with the 
non-l inear predictions obtained from CAMB (|Lewis et al.l 
2000) up to k ~ 2.0/iMpc -1 , which sets the resolution 
scale of our power spectrum measurements. As expected 
from WNLWF, the Gaussian component of WNLWF pre- 
serves the power on linear scales - up to k ~ O.lMVIpc -1 
- while signals from trans-linear and non-linear scales are 
mostly transferred to the non-Gaussian contribution of the 
WNLWF decomposition, which explains the drop in power. 
At the same time, we note that the power spectrum after 
WNLWF actually traces quite well the linear predictions, to 
within a factor of two, at all scales. 

The density reconstruction algorithm also has a signifi- 
cant impact on the shape of the power spectrum, as particles 
are pumped out of the gravitational potential. As a result, 
small scale power is directly transferred to larger scales, as 
seen in the figure. Interestingly, the turning point between 
these two effects is also locate close to k ~ 0.1/iMpc - . 
This common feature to both Gaussianization technique is 
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Figure 1. Projections through a thickness of 50 cells of one of the realizations. In each panel, the side is 322.36/i _1 Mpc, and the image 
contains 1024 2 pixels. Top left is the original field, top right is the field after linear density reconstruction, middle left is the wavelet 
filtered field (Gaussian part), middle right is the result of wavelet filtering the reconstructed field. The non-Gaussian part of the wavelet 
filtered fields are shown in the bottom panels with (right) and without (left) density reconstruction. To ease the visual comparison, each 
panel shows the same ovcrdcnsity range and saturates for denser regions, i.e. all pixels with 8 > 30 are black. 
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Figure 3. Cross power spectra of a single density field, con- 
structed from a random separation of the haloes onto two dis- 
tinct fields, whose Fourier transform are combined. The dashed 
line is the halo power spectrum of the full population, the dotted 
line is the power spectrum of the particles, the straight line is 
the Poisson noise estimate, the open symbols represent the cross 
spectrum of the two randomly selected populations, and the thick 
solid line is the shot noise subtracted power. 



explained by the fact that both attempt to pump back infor- 
mation from scales in the trans-linear regime, whose lower- 
k limit indeed corresponds to this mode. When looking at 
the halo measurements, we observe that the original and re- 
constructed power spectra are dominated by shot noise at 
scales smaller than k ~ 1.0/iMpc -1 . This noise, however, is 
strongly suppressed by the wavelet filter. 

In practice, a common way to deal with Poisson noise 
is to compute the cross power spectrum between two popu- 
lations randomly selected out of the original catalogue. The 
shot noise is eliminated in this operation, and the signals left 
behind are stronger on large and intermediate scales. Small 
scales are typically anti-correlated due to the 'halo exclusion' 
effect, which is a result of our halo-finder that collapses all 
the structure to a single point, and leaves the surrounding 
region empty. Although precise and robust, this procedure 
is hard to apply to all the cases under study in this paper, 
since the density reconstruction algorithm looses efficiency 
if we under-resolve the small scale structures. Instead, we 
use another common approach which consists in subtract- 
ing from the measured power spectrum an estimate of the 
shot noise, defined as P a hot = V 3 /jV ha i os . To illustrate this, 
we show in Fig. [3] a comparison between halo power spec- 
trum from the full catalogue, the cross spectrum, and the 
shot-noise subtracted power. We observe that the two shot 
noise subtraction techniques agree up to k — 0.3/iMpc -1 , 
beyond which the P(k) — P s hot(k) approach looses power in 
comparison; by k = 1.0/iMpc - , it is smaller by a factor of 
2.4. This difference between the two shot noise subtraction 
techniques means that our approach is not optimal, and that 
the results on the smallest scales are not as robust as one 
would wish. 



3.2 Covariance matrices 

The two Gaussianization techniques that are discussed in 
this paper both attempt to bring cosmological information, 
or degrees of freedom, back to the power spectrum. Conse- 
quently, the covariance matrices of the Gaussianized fields 
will be more diagonal. The top left panel of Fig. [3] shows 
the cross-correlation coefficient matrix of the original parti- 
cle fields in the upper triangle, and the wavelet filtered ones 
on the lower triangle. To ease the comparison between the 
figures, we show in the main figure of each panel the posi- 
tive components only, and present in the insets the negative 
entries. There is a mild anti-correlation (less than 10 per 
cent) in some matrix elements of the original fields, which 
comes from residual noise in the largest scale. This is a mild 
effect that has very little impact, hence we do not attempt 
to correct for it. The top right panel shows, on the upper 
triangle, the results after a density reconstruction has been 
applied, then, on the lower triangle, the measurements after 
both technique have been executed. The off-diagonal ele- 
ments of the covariance matrix are reduced by 20-30 per 
cent by both Gaussianization techniques. We see that those 
two techniques combine well and reduce to a minimal value 
the correlation between the modes. 

The bottom panels of Fig. U show the same mea- 
surements, when carried on halo fields. The wavelet fil- 
ter produces a band of negative elements, correlating the 
k > 1.0/iMpc -1 - Poisson dominated - with all scales. This 
anti-correlation does not carry any physics about the signal 
in it, hence these scales should be left out or carefully in- 
terpreted in future analyses. Most off-diagonal elements are 
about 30 per cent less correlated than in the unfiltered ma- 
trix, showing that wavelet filtering is also very efficient on 
halo fields. We also observe that the density reconstruction 
algorithm has very little impact on the correlation of the 
halo measurements. This is caused by the fact that haloes 
are non-overlapping by construction, hence the region of 
exclusion prevents an accurate construction of the gravita- 
tional potential This is a strong limit of the technique, as 
real galaxy data behave much more like halos than parti- 
cles. However, higher multiplicity in the galaxy population 
of many haloes is likely to improve the construction of the 
gravitational potential, hence we can expect a modest gain. 
When we combined both technique, most of the diagonal- 
ization comes from wavelet filter. This can be seen visually 
by comparing the lower triangles of both bottom panels. 

3.3 Fisher information 

When we extract the Fisher information from the covari- 
ance matrices presented above, we expect the original par- 
ticle fields to exhibit th e global shape first measured in 
iRimes fc Hamiltonl (|2005l ). Namely, the information should 
follow the Gaussian predictions on large scales, then reach 

8 We could of course improve the performance of the reconstruc- 
tion technique by using the gravitational potential measured from 
simulated particles, which we have at hand. Even in a data set, it 
is in principle possible to combine independent measurements of 
the potential, obtained say with weak lensing tomography. This 
is an interesting avenue that is, however, beyond the scope of this 
paper. 
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Figure 2. (top:) Power spectra of the original and Gaussianized fields, from simulated particles (symbols + solid lines) and haloes 
(symbols + dashed line, in red on the on-line version). The linear and non-linear predictions from CAMB are shown by the thick solid 
lines, and the Poisson noise corresponding to the halo population is shown with the thin dotted line. We observe that the wavelet filtered 
particles densities trace the linear CAMB predictions, at most within a factor of two. (bottom:) Fractional error between the curves of 
the top panel and the non-linear prediction from CAMB. We observe that the particle power spectrum deviates by more than 10 per 
cent for k > 2.0feMpc — 1 , which sets the resolution limit of our simulations. This scale is represented by the vertical line in both panels. 
We observe that in linear regime, the wavelet filter preserves the agreement with the predictions, whereas the reconstruction tends to 
increase the power spectrum by about 20 per cent. This is not a surprise since one effect of the algorithm is to transfer power from 
small to large scales. In the non-linear regime, however, the power spectrum is highly suppressed by the wavelet filtering process, which 
factorizes the structures into the non-Gaussian contribution. We measure a linear bias of about 1.2 in all halo measurements. The original 
and reconstructed halo power spectra are shot noise dominated at scales smaller than k ~ 1.0/iMpc - , a scale that is strongly suppressed 
by the wavelet filter. 



a trans-linear plateau where the gain is very mild as one 
increases the resolution of the survey, then hit a second rise 
on scales smaller than about 1.0/iMpc -1 . We first see from 
Fig. [5] that we are a ble to recover those results, plus those 
of Ngan et al.l ([2011) , which showed that the density recon- 
struction algorithm can raise the height of the trans-linear 
plate au by a factor of a few. We also recover the results 
from I Zhang et all (|201lT) and obtain a similar gain with the 
wavelet non-linear Wiener filtering technique. 

As mention ed in the introduction, it was shown by 
lYu et all l|2012f ) that different Gaussianization techniques 
do not always combine well. In the current case, however, 
we observe that on all scales, the Fisher information from 
the combined techniques are larger than the sum of the two 



separate contributions. For k > O.B/iMpc" 1 , notably, we are 
able to extract more than ten times the Fisher information 
of the original particles fields, whereas individual techniques 
offer a recovery of about a factor of four. This symbiosis ef- 
fect grows larger as one goes to smaller scales. 

When considering the halo fields, we observe in Fig. [5] 
that the density reconstruction technique, taken alone, has 
little impact on the recovery of information, due to a poor 
modeling of the gravitational potential. In contrast, wavelet 
filtering recovers five times more information by the time we 
have reached k = 1.0/iMpc - , before shot noise subtraction. 
The Poisson noise is a non-Gaussian effect, which also sat- 
urates the Fisher information. A hard limit one can think 
of is the following: the number of degrees of freedom can 
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Figure 4. (top-left) Cross-correlation coefficient matrix associated with the particle power spectra. The top triangle represents mea- 
surements from the original matter fields, while the lower triangle elements are from fields after the NLWF has been applied. The inset 
quantifies the amount of anti-correlation between the measurements, (top-right) The top triangle represents measurements from the 
reconstructed matter fields, while the lower triangle are from fields that are first reconstructed, then wavelet filtered - still using parti- 
cles as tracers, (bottom-left:) Cross-correlation coefficient matrix associated with the power spectrum measurements from the simulated 
haloes. The top triangle represents measurements from the original halo fields, while the lower triangle are from fields that are wavelet 
filtered, (bottom-right:) The top triangle represents measurements from the reconstructed halo fields, while the lower triangle are from 
fields that are first reconstructed, then wavelet filtered. 



not exceed the number of objects in our fields of view. We 
therefore plot the (non-Gaussian) Poisson noise limit as a 
flat line corresponding to the halo number density, and ob- 
serve that the original halo Fisher information approach but 
never exceed that limit. Wavelet fields, however, reduces the 
Poisson noise significantly, hence allows the information to 
reach higher values. We also see that shot noise subtracted 
Fisher information curves show a lower information recov- 
ery, which means that the number density needs to be high 
enough in order to maximize the recovery. As mentioned in 
section [2~T1 the halo density is about eight times larger than 
current spectroscopic surveys. Next generation experiments 
and current photometric redshift surveys have a much larger 
number counts, hence the corresponding Poisson noise limit 
will be much higher. 



4 DISCUSSION AND CONCLUSION 

This paper explores the recovery of Fisher information with 
the combined use of two Gaussianization techniques, and 
show that wavelet non-linear Wiener filtering and density 
reconstruction can extract an order of magnitude more in- 
formation than the original fields. We also reproduce the 
calculations on halo catalogues and find that 1) the den- 
sity reconstruction has only a mild impact on its own, 
due to a poor modelling of the gravitational potential, 2) 
wavelet filter recovers about five times more information by 
k — l.OftMpc , and 3) the combined techniques recovers 
about three times more Fisher information than in the orig- 
inal fields by k > 0.7/iMpc -1 , even after shot noise sub- 
traction. Interestingly, we find that in both matter tracers, 
the recovery of the combination of the two Gaussianization 
techniques is larger than the sum of the individual contri- 
butions. 
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Figure 5. (top :) Cumulative information contained in the dark matter power spectra of the original and Gaussianized fields from particles. 
As in Fig. [2] the dots represent the original fields, the open circles show the results after our density reconstruction algorithm, the squares 
correspond to the wavelet filtered fields, and the stars represent a combination of both techniques. The analytical Gaussian (i.e. linear) 
Fisher information curve is shown with the thick solid line. These two Gaussianization techniques are shown to work in conjunction, such 
that on all scales, their combined effect recovers the largest amount of information. For k > 0.6/iMpc - 1 , the improvement on particles 
is more than an order of magnitude, (bottom :) Ratio of the lines presented in the top panel with the original fields. 
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